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LONG-TERM STORAGE CAPACITY 
RESERVOIRS 


solution the problem determining the reservoir storage required 
given stream, guarantee given draft, presented this paper. For ex- 
ample, long-time record annual total discharges from the stream avail- 
able, the storage required yield the average flow, each year, obtained 
computing the cumulative sums the departures the annual totals from 
the mean annual total discharge. The range from the maximum the mini- 
mum these cumulative totals taken the required storage. 

The results the investigations given here were applied computing the 
storage required the Great Lakes the Nile River Basin, the length time 
chosen basis the estimates storage being 100 years. These reservoirs, 
combination with the other projects described, would enable the Nile 
developed for irrigation the fullest possible extent. relation has 
been derived between the storage capacity the main stream and the amounts 
required the tributaries produce its equivalent. 

thought that the general theory may have other applications than the 
design reservoirs for the storage water. 


INTRODUCTION 


number projects has been investigated which will enable Egypt and 
the Sudan develop irrigation from the Nile River its fullest 
The broad plan involves storage water from good years for use bad ones, 
and necessitates reservoirs sufficient capacity meet the shortages that 
might occur during century. Much research had been done the capacity 


comments are invited for publication; the last discussion should submitted 
September 1950. 

Consultant the Ministry Public Works, Egypt; formerly Director-Gen., Physical 
Dept., Cairo, Egypt. 

Nile Vol. VII, Hurst, Black, and Simaika, Physical Department 
Paper No. Ministry Public Works, Cairo Govt. Press, Cairo, Egypt, 1946. 
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such reservoirs and the results were used determine the size the reservoirs 
required the Great Lakes the Nile Basin for what was called 

Some preliminary results investigations were released which 
were amplified 1946;? but neither case was there adequate statement 
the theoretical investigation. The latter has been reserved for this paper. 

Storage problems were treated first the late Allen ASCE, 
1914, when studied the discharge thirteen American rivers, dealing with 
them mainly graphical manner, probability paper. The longest record 
was one years, but longer sets observations were made combining 
records from several rivers. The objection this procedure and the uncer- 
tainty extending the normal Gaussian frequency distribution the future 
were recognized Mr. Hazen, but thought gave the best results then 
obtainable. 

The work was extended the late Charles ASCE, similar 
graphical methods. Much Mr. Sudler’s work was based artificial records 
which information taken from short record stream was extended 
period 1,000 years writing, say, annual runoff values cards and 
then shuffling and drawing card from these 1,000 times. Again the actual 
data were necessarily confined short periods since stream measurements the 
United States had rarely been made for long years. Thus, the results 
are mainly based few short-period records natural river flow extended 
long periods random repetition. clear that the process can give very 
little information about what may happen over long periods and may com- 
pletely misleading since excludes values higher lower than those which 
have already occurred. This was recognized some those who contributed 
the discussion Mr. Sudler’s paper. Moreover, the graphical methods 
employed not lead any concise and understandable presentation 
results. 

the present investigation the methods the theory probability are 
applied, but the data are long-period records natural events which are shown 
have certain similarities, although they may records river discharge, 
rainfall, temperature, annual growth rings trees, annual deposits clay 
lakes. The average results can expressed quite simply two equations 
which apply the storage other substances well water. shown 
that, regard storage, records covering such short periods years 
years may very misleading, but their usefulness can extended 
analyzing large number phenomena which have been recorded for periods 

The problem what storage required stream give certain 
minimum discharge was first investigated the case Lake Albert (in the 
Belgian Congo and the Uganda Protectorate) determine the size over- 
year storage reservoir which would required order equalize the outflow 


Nile Basin,” Vol. Hurst and Phillips, Cairo Govt. Press, Cairo, Egypt, 1938, 86. 
Provided Impounding Reservoirs for Municipal Water Allen Hazen, 
Transactions, ASCE, Vol. LXXVII, December, 1914, 1539. 
Required for the Regulation Stream Charles Sudler, ibid., Vol. 91, 1927, 
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over number years. leads statistic which the range from maxi- 
mum minimum the curve obtained plotting the cumulative totals 
the departures from the mean the annual discharges taken order. Such 
curve cumulative totals sometimes called mass curve. Values 
have been computed for long series observations rainfall, temperature, 
and other natural phenomena, and these have been filed for reference the 
Engineering Societies 

letter symbols introduced this paper are defined 
where they first appear, the text, illustration, and are arranged alpha- 
betically the Appendix, for convenience reference. few specialized 
terms used the paper are defined follows: 


Century computed storage capacity reservoir, re- 
quired the Great Lakes the Nile Basin, necessary meet the shortages 
that may occur 100 years (Section 1). 

Continued Sum.—Successive cumulative totals (Section 1). 

Cubic Meter.—English units, 0.0008107 acre-ft. 

Equatorial Lakes.—Lake Victoria (in Uganda, Tanganyika, and Kenya), 
Lake Albert (in the Belgian Congo and Uganda), and other bodies impounded 
water the southern Nile Basin, the equator; also referred the 
African simply the (Section 2). 

Gaussian Frequency normal frequency curve first advo- 
cated the German mathematician, Karl Freidrich Gauss, who maintained 
that all skew distributions tended become normal the number events 
was increased and that skewness was evidence inadequate data. 

Meter.—English units, 3.28083 ft. 

thousand million things—that is, 

Varves.—Annual layers silt deposited lake other body still 
water. individual layers differ thickness and character, succession 
such layers forms characteristic group which .can identified con- 
temporaneous deposition whatever deposit may found. thus 
possible, combining different sections, measure the time involved the 
deposition the entire group sediments and construct time scale 
manner similar that employed the study annual rings trees (Webster’s 
New International Dictionary the English Language). 


STATEMENT THE PROBLEM 


Given reservoir large capacity, how can its outflow regulated eco- 
nomically meet the deficiencies low years, and what the capacity required 
guarantee certain minimum discharge these years? These are the 
questions. 

concrete example illustrate the methods, consider the case Lake 
Albert. the past annual discharges from the lake are given, the size 
reservoir required maintain the maximum possible steady discharge during 
the period can determined. This the mean discharge for the period, and 
the capacity the reservoir obtained adding, algebraically, the de- 

39th Street, New York 18, 
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partures from the mean form series accumulated totals. The difference 
between the highest and the lowest these accumulated totals—that is, the 
range R—is (a) the maximum accumulated storage when there never 
deficit, (6) the maximum accumulated deficit when there never any storage, 
(c) their sum when there both storage and deficit. This range the 
storage required maintain the average discharge. 

Increased losses due storage are disregarded because, unless they are 
small, the site not suitable for over-year storage. reservoir like the 


Annual out flow. 
Continued sum 

AA, Mean annual out flow. 


metres. 


1903 1910 1920 1930 1940 1946 


Aswan Reservoir Egypt there heavy evaporation uncompensated any 
rainfall, and the surface area the reservoir several times that the natural 
river. Therefore, there appreciable extra loss result storage, 
that, even there were the capacity Aswan, would unsuitable for over- 
year storage. the case the Equatorial Lakes their size makes possible 
store great volume water with only small proportional increase 
surface area. Moreover, evaporation and rainfall are nearly equal; in- 
creased evaporation loss due increased surface compensated increased 
gains from rainfall. 


a 
-40 
-60 
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Fig. shows the annual total discharge from Lake Albert based recorded 


lake levels and gage-discharge curves. 


The curve accumulated departures 


from base near the mean also shown. The actual values and the method 
computation are shown reference Table 


The discharges, Col. are 
usually given one tenth 
milliard cubic meters but the last 
digit not certain. For this com- 
putation the discharges are written 
the nearest milliard and this re- 
duces most the departures, 
single digit, which enough for 
purposes computation. The omis- 
sion the last digit makes sig- 
nificant difference the mean the 
standard deviation, but lessens 
considerably the labor computa- 
tion. The computation arranged 

The mean the series first 
computed, and then base chosen 
near the mean. this case the 
selected base from which 
compute departures, being the near- 
est whole number the mean, and 
keeping all the departures one 
two digits. 

After writing the departures the 
positive and negative values 
(Cols. and are added separately. 
this case the sum the positive 
122 and the negative, 114. The 
difference and 8/N should 
the difference between the mean and 
the base. difference gives 
check the correctness the de- 
partures, and also their continued 
sum 2d, (Col. 6). the departures 


were reckoned from the mean 


final algebraic sum would zero; 
otherwise, (mean minus 
base), (24.19 24), and, 


course, this also the algebraic sum the totals the positive and negative 
departures, thus affording check the continued sums. 
series observations advisable compute the totals the columns, 


TABLE From 
ALBERT (BELGIAN 


PROTECTORATE) MILLIARD 
METERS 810,700 


(1) (2) (3) (4) (5) (6) 
1943 


1,040 +122 —114 2,310 


each page check the work piece piece. 


With long 
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computing the standard deviation, the mean square deviation must 
0.04; and 7.33. This the only part the work which 
not self-checking and must checked repetition. The last digits the 
square are omitted, which means this case that the squares are given the 
nearest ten which near enough provided there are more observations. 
The effect the approximation about this instance. 

The departures are reckoned from base which not the mean, and, the 
continued sums these are plotted, simple graphical method will yield the 
continued sum departures from the mean any other base. the present 
case the mean 24.2 and the base, 24. departure from the base 
therefore 0.2 greater than the corresponding departure from the mean, that 
the sum the first departures 0.2 greater than the sum the first 
departures from the mean. The sum all the departures from the base 
0.2. Thus, the continued sum any number departures from the 
mean can found joining the point representing the last continued sum 
the origin and taking this line new axis. This useful construction for 
finding the continued sums from any base, when the sums referring some 
base have been plotted. 

Fig. the continued sums departures from base are plotted and 
the new axis giving continued sums departures from the mean drawn. 
The range, the summation curve departures from the mean 
milliards. 

Examples other uses the graphic method obtaining continued sums 
are given, subsequently, Figs. and 

the case Lake Albert the accumulated sum departures rose 
maximum 1919 (point Fig. and then fell zero 1946 (point X), 
that the next period would begin with stored water. Had the operation 
been begun 1920 there would have been deficit start with. This fact 
shows the necessity making reservoir some time before will required 
that can partly wholly filled before its storage needed. clear 
from Fig. that the storage required maintain the mean discharge over the 
period the range from maximum minimum the accumulated departures. 
also clear that, the capacity the reservoir make 
the maximum accumulated deficit, when this has been done, the reservoir will 
empty and there will margin stored water for the future. Moreover, 
the engineer regulating the flow after the event, with full knowledge what 
has happened. The problem the required capacity not simple one and 
cannot solved adding departures from the mean, considering only Lake 
Albert, but must dealt with wider basis. The discharge the main 
Nile River Aswan may taken further illustration the difficulties. 

Table for example, based three widely differing assumptions: First, 
that computation was made 1909, using years record; second, that 
computation was made 1946, using the succeeding years record; and, 
third, that computation was made 1946, using the total period years 
record, The difference between the results these computations shows 


4 
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how misleading the simple process can be, both the size the reservoir 
needed and the benefit expected fromit. Thus, clear that more 
extended investigation required, and this can only made the examina- 
tion many sets river discharges extending over long periods, other 
data with similar characteristics. the same time some promise solution 
possible exploring theory. 


TABLE 2.—Comparison COMPUTATIONS FOR 
THE Main Nite 


Description 1871-1908 1909-1945 1871-1945 
Average annual discharge, milliards 103 
Range accumulaied departures, milliards 201 476 


The first part this paper attempts find theoretical relation among the 
range, the summation curve departures from the arithmetic mean, the 
variability the phenomenon expressed its standard deviation, and the 
number years observation, the assumption that the observations 
are random events following the normal Gaussian frequency distribution. This 
course was adopted because many natural phenomena have frequency dis- 
tributions which are approximately normal the order occurrence dis- 
regarded, and was thought that the examination this case, which first 
approximation those occurring nature, would indicate the form relation 
sought. Gauss’ normal frequency curve (which also the curve for the 
frequency distribution errors observation) gives the frequency depar- 
tures quantity from its arithmetic mean relation their magnitudes. 
number values the quantity the equation the curve 


Incidentally, the assumption that the discharge river random event 
the basis the works Messrs. and under this same assump- 
tion theoretical solution has been found. 


FREQUENCY DISTRIBUTIONS RANDOM EVENTS 


well known that the Gaussian normal frequency curve (Kq. fits 
many distributions events which success failure are equally likely and 
which each event independent the next—for example, the frequency 
accidental errors set observations. also well known that approxi- 
mation this curve given the terms binomial expansion. When 
set coins tossed, the probability the occurrence heads and r)- 
tails any throw mC, which the number combinations 

the set tossed times (the number being fairly large), the average 
frequency occurrence heads and tails, head and 1)-tails, etc., 
given the terms the binomial expansion: 


| 
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increases, Eq. approximates the normal curve, and the ap- 
proximation fairly good; fact, for many purposes good enough. 

Table shows how the number heads would expected distributed 
the average when set coins tossed 1,024 times—that is, the binomial 
distribution. The distribution derived from the normal curve, and one actually 
obtained trial are also given. 


TABLE FREQUENCIES 


NuMBER 


Distribution 


the normal distribution the numerals represent the nearest integer. The 
relation between the binomial and the normal distribution has been treated 
Udny who defines the normal distribution the limit the 
binomial distribution when indefinitely increased. For theoretical pur- 
poses the computer can use whichever form the more convenient. 


CALCULATION THE RANGE CONTINUED SuMS DEPARTURES 
STORAGE FOR FREQUENCY DISTRIBUTION 


Assume number consecutive values element whose variation from 
its mean distributed normally, the binomial expansion already dis- 
cussed, but where there correlation between successive values groups 
values, tossing coins. For example, consider the definite case tossing 
purposes each toss the coins recorded the number heads minus 
the number tails, and the departures this quantity from its mean are 
added find the range between the maximum and the minimum this 
summation curve. The range thus found, from tosses set coins, 
may differ slightly from the range found the summation heads minus 
tails from tosses single coin. This can easily seen tossing 
set coins one after the other, and recording individual results. Then the 
sums heads minus tails are the same the end every set ten whether 
individual tosses are considered not. the continued sum 
plotted using single tosses, every tenth point being marked. The curve ob- 
tained taking the tenth points would close approximation the curve 
individual points, but the latter would usually have slightly greater range. 
For example, suppose that the maximum the curve was produced set 
heads and tails which would put the end-point the set two units above 
the end-point the preceding set. the heads occurred and then the 
tails, there would point the curve recording individual sums heads 
minus tails which would six above the end-point the previous set—that is, 


Introduction the Theory Udny Yule, Charles Griffin and Co., London, 
Ed., 1912, 317. 
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the range would increased four. For the present this difference 
ignored, but will considered later. 
the case selected for analysis coins are tossed times, heads and 
tails gains and losses are equal the end the trial, and written for 
The number orders arrangements which gains and losses can 


occur 


time. Among these orders there are instances which, some 
point the process tossing, losses exceed gains 

There are three cases: 


Case Condition No. arrangements 
Losses never exceed gains 
(a) 
Gains never exceed losses 
(c) Losses never exceed gains 
and gains never exceed 
n+1 
losses 


case (a) the range the maximum value gains minus losses; case 
the range the maximum value losses minus gains; and case (c) the 
range the sum the maximum gains minus losses, and the maximum 
losses minus gains. 

Summing all values the range obtain mean value, and writing for 
number gains and for number losses, 


onC'n 


(4) 


find the mean range see Table There are arrangements 
which, some point, gains exceed losses one more; and there are 


TABLE DETERMINE THE MEAN RANGE 


Gains minus losses No. arrangements Products 
1 wnCn41 = onCn42 1 (onCn41 inCn+2) 
2 onCn42 onCn4s 2 (onCn42 onCn43) 
3 — 3 (onCn42 — 


n n (2nC an) 


Sum onC'n41 onCn41 + anCn42 + + onCon 


“Choice and Chance,” Whitworth, Cambridge, England, 4th Ed., 1886, Chapter Priority. 
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which gains exceed losses two more. there are 
arrangements which gains exceed losses one. 

Referring to Table 4, 1 + + onC'2 + + + + 
The arrangements Table include all except those which gains never 


exceed losses, whose number (or and which the 


maximum gains minus losses Hence, the mean maximum value 
2nvn 
minus gains. Therefore, the mean range 


gains minus losses and this also the mean maximum losses 


Since large, Eq. can simplified the approximation for James 
Stirling which follows: 


Hence, 
and 


since the product large. 

Thus the average range the accumulated sums (number heads minus 
number tails), when coins are tossed times being large), in- 
creases the square root the number tosses, exactly like the accumulated 
error line leveling. 

The standard deviation the binomial distribution produced tossing 
coins 


This the standard deviation the number heads (r) tails 
The standard deviation the number heads minus the number tails 


coin, and little larger than the range obtained from tosses set 
coins. The average difference between the two computations for range 
can found follows: 

Toss coin large number times and plot the continued sums heads 
minus tails, marking given number them off sets 2m. The set 


RESERVOIR CAPACITY 


that contains the maximum sum selected for scrutiny. there are 
heads and tails the end-point the set higher than the end-point 
the previous set. The number orders which heads and tails can occur 

The number orders which some stage the tossing the number 
heads shall exceed the number tails more the end 
the set, heads exceed tails The number orders 
which some stage heads exceed tails more than 
all these cases the maximum the summation curve individual tosses 
will exceed the maximum the curve derived from sets tosses. The 
mean value the excess can obtained procedure similar that 
already employed finding the mean range the summation curve. The 
result this 


m)! 


There similar effect the minimum the summation curve, that 
the average range the summation curve individual tosses greater than 
the range the curve from sets tosses 2d. Some values are given 
increases. 


VALUES 


Magnitude 
. a 
Table 6(a) gives the results 4.1 4.3 4.5 4.8 


actual experiment which sixpences 
were shaken box and thrown 
onto table. They have been ana- 
lyzed some detail give practical example the foregoing theory. 
Their frequency distribution has been given Table which there were 
1,024 throws coins which resulted distribution close the 
theoretical. 

The observations were recorded they occurred and the continued sums 
were computed from 1,000 observations. this case (because, the 
finish, heads and tails are nearly equal) the actual values (heads minus tails) 
are taken departures, and not the departures from the mean (Col. Table 
6(a)). the finish there are more heads than tails; the mean value 
per toss 0.074. This corrected graphically the aforementioned 
method. Fig. discussed subsequently, shows the result one the 
experiments. 

each case (Cols. Table corrected what would have been 
the departure had been reckoned from the mean. The theoretical value 
3.16 for binomial distribution. 
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Table 6(a) shows that tends constant and the values from 
1,000 trials (Cols. agree closely with the theoretical values (Cols. 5). The 
trials illustrate point not examined the theoretical work, and that the 
variability and from set set 100 throws. 

The maximum variations, percentage departures from the mean, are 
follows: 11%; 50%; and 39%. The standard deviation (or 
standard error) for set 100 throws about 7%. 

Another device for obtaining normal frequency distribution what may 
called probability pack cards (Table this pack the cards are 
numbered +1, —1, +3, —3, +5, —5, +7, +9, and —9, and the numbers 


TABLE witH RANDOM EVENTS 


DEFINITIONS COLUMN HEADINGS 
Primary Observations: 
Table Excess heads minus tails —t. 
Table Sum recorded numbers 
Table 6(c), Excess even numbers over odd numbers 
Col. Standard corresponding primary observations (Col. 
Col. Range accumulated deviations 


from observed data. 


Col. Computation 


Col. Theoretical values data Cols. corrected for the effect tossing sets. 


1,000 Times 1,000 
rials 


each kind are proportional the corresponding ordinates normal fre- 
quency curve. They are: one’s; three’s; and and respectively, 
the others. The approximation these numbers fairly close. The cards 
are first well shuffled and then cut, and the number the exposed card 
recorded. The cards are reshuffled slightly and cut again, and on. The 
numbers recorded may taken corresponding observations quantity 
whose frequency distribution conforms the normal Gaussian curve. 

This quicker process than tossing and counting coins. toss coins 
100 times required about min, whereas shuffling and cutting for 100 cards 
required min. Table records the results trial cutting 1,000 times. 
Col. Table 6(a), each case corrected what would have been 
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the departures had been reckoned from the again, 
Cols. approximates the value obtained theoretically for the case when 
large (Col. 5). 

third set chance events was examined (Table 6(c)). From time 
time newspapers announce the numbers bonds that have been drawn for 
redemption. collection unconnected random events was made dividing 
the numbers bonds into groups ten, follows: 


(b) 


— 


Number Even Digits Minus Number Odd Digits 


The left-hand column group was ignored only few different digits 
are likely occur it. The number even digits less the number odd 
digits columns (a) (c), (d) (e) (f), are recorded. This difference 
yields sequence numbers exactly analogous the sequence previously 
obtained tossing coins, and the labor producing much less. 

this method sequence 1,000 random numbers composed the even 
numbers from +10 was made and this was used before. Table 6(c) 
gives the results this trial. 

each case (Col. has been corrected what would have been 
the departures had been reckoned from the mean. The theoretical value 
(Col. has been corrected for the effect tossing sets. 

Although the mean value found from the three sets random 
events Table very close the value found mathematical analysis 
(Cols. compared Cols. yet individual values for sets 100 observations 
vary considerably. The extreme ranges the three cases are follows: 


Table 


The mean value from the thirty sets 100 observations 1.27 
and the standard deviation 0.32. The relation between and for 
random events having been established theoretically and having been found 
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agree with practice, now necessary find what respects natural 
phenomena are similar random events and what respects they differ. 


CHARACTERISTICS NATURAL PHENOMENA 


The phenomena that chiefly concern the subject this paper are river 
discharges. Unfortunately records these rarely cover periods more than 
years, although levels (which are usually indication the discharge) have 
been recorded over longer periods. The discharge river closely related 
the rainfall its basin, and rainfall records for many places cover more than 
100 years. The investigation began with these phenomena. 

The main similarity river discharges, levels and rainfall, and random 
events that their frequency distributions conform approximately the 
normal curve, although they are usually slightly skew. rule the observa- 
tions are too few give very regular 
distributions, but the flood levels 
values. (nearly 1,000 them), the Nile 
River Cairo, Egypt, recorded 
the Roda gage, give very regular 
distribution which well fitted 
normal curve. This distribution 

Although random events 
groups high low values 
occur, their tendency occur 
natural events greater. This 


per thousand. 


25 -20 -10 -O5 0.5 20 «(2.5 
from mean metres. 


Between 641 A.D. 
1946 A.D. 


the main difference between natural 
and random events, example 
which the discharge the Nile 
River Aswan. The long series 


records flood levels Cairo shows the same There 
obvious periodicity, but there are long stretches when the floods are generally 
high, and others when they are generally low. These stretches occur without 
any regularity either their time occurrence duration. 

the investigations Messrs. Hazen and Sudler the frequency character- 
istics river discharges are assumed like those random events, and 
this the common assumption estimating probability high low floods. 
will shown that this assumption must used with caution and, cases 
which storage over long periods concerned, only approximation. 


FoR PHENOMENA 


The investigation began experimentally, before any attempts find 
mathematical theory were made, finding for any available long series 
river discharges. Such series were scarce and the work was extended 
rainfall data, which are more plentiful. When was found that rainfall data 


Records the River Jarvis, Transactions, ASCE, Vol. 101, 1936, 
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gave results similar the river data the work was extended other long series 
natural events. 

For example, Fig. shows summation curve departures for the maxi- 
mum levels lake the Dalalven River Sweden. There are 176 years 
observations, divided into periods years, years, and 132 years; and 
the ranges, the summation curves for the departures from the mean 
each these periods are found. 

The sums departures from the base are plotted. For the sum 
departures from the mean for the entire period the curve referred the 


1000 


500 


1764 7800 1850 7900 1940 


LAKE THE DALALVEN RIVER SWEDEN 


axis OD. Similarly, the summation for the first years referred line OA, 
for the second years line AB, for the first years line OB, and on. 
The range, the last period PM. 


The theoretical investigation shows that proportional the case 
random events, and this must approximately the case for any events that 
have normal nearly normal frequency distributions. The device express- 
ing deviations from the mean terms common one, derived from the 
expression for the normal curve; and, since function these deviations, 
has been calculated. Complete tabular records the computation in- 


cluding mean values and are the complete original 
manuscript, cataloged for reference the Engineering Societies 
Seventy-five different phenomena were used, and 690 values were com- 
puted for different values 
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summary the results given Table which the phenomena are 
divided into groups according their nature and the value and log 
plotted against log Fig. 

general, has not been computed for periods less than years be- 
cause, for very short periods, has little practical value although may have 
some formal value. For has the exact value, unity. The longest 
period used 2,000 years and this for clay layers Lake Saki the Crimea. 
Fig. which the means small groups observations the same kind 
and length period have been plotted, shows that log linear func- 


7.5 1.54 0.85 0.68 
River levels, discharges, and runoff....... 8.9 1.65 0.94 0.70 
13.1 1.79 1.08 0.72 
108 16.4 2.02 1.19 0.69 
105 19.6 2.02 1.27 0.75 
208 36.5 2.32 1.54 0.77 
309 53.9 2.49 1.72 0.79 
420 60.3 2.56 1.77 0.78 
613 2.79 1.89 0.77 
716 104 2.85 2.01 0.79 
820 122 2.91 2.08 0.80 
927 129 2.96 2.11 0.79 
1,040 130 3.02 2.12 0.78 
Rainfall stations with three groups 8.2 1.57 0.91 0.72 
121 22.2 2.08 1.31 0.74 
Temperature and pressure stations with 7.4 1.55 0.85 0.68 
110 18.0 2.04 0.71 
8.8 1.63 0.93 0.70 
Temperature stations with values R.. 15.1 1.94 1.17 0.71 
132 21.0 2.12 1.32 0.72 
175 27.1 2.24 1.43 0.74 
13.6 1.69 1.11 0.80 
100 26.6 2.00 1.43 0.84 
Annual growth tree 200 45.0 2.30 1.64 0.82 
300 79.9 2.48 1.90 0.87 
462 76.6 2.66 1.87 0.79 
900 187.0 2.95 2.27 0.86 
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7.—(Continued) 


No. 
Phenomena years Log Log— 
10.9 1.70 1.02 0.73 
Thickness annual layers mud; Tamis- 100 21.3 2.00 1.31 0.77 
kaming, Ont., Canada, and Moen, the 200 42.7 2.30 1.58 0.79 
Sogne District, Norway............... 300 82.5 2.48 1.90 
400 126 2.60 2.09 0.91 
550 115 2.74 2.00 0.82 
1,100 181 3.04 2.19 0.80 
1.98 1.30 0.77 
9.7 1.70 0.98 0.70 
Thickness annual layers mud, Lake 100 15.3 2.00 1.17 0.69 
500 47.9 2.70 1.66 0.69 
1,000 84.0 3.00 1.91 0.71 
2,000 179.0 3.30 2.24 0.75 
12.4 1.77 .06 0.72 
Sunspot numbers and wheat prices....... 124 22.1 2.09 1.34 0.75 
237 16.9 2.36 0.60 


Rainfall stations with one value includes temperature one station. Rainfall stations with 
two groups values includes temperature and one pressure. ranges from 120. 


tion log and can expressed 


each main group observations line drawn through the point log =0, 
log 0.30 2), and the center gravity all the observations 
good fit. The slopes these lines (K) vary from 0.69 0.80. 

Whether there any theoretical significance the fact that approxi- 
mately not known. The value has been calculated for each the 
690 individual values and the frequency curve shown Fig. has been 
drawn. The data for this curve are follows: 


Description 
Number values................ 690 
Range individual values................ 0.46 0.96 
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BETWEEN RANGE SUMMATION THE STANDARD DEVIATION 
AND THE YEARS 


Fig. shows that there slight skewness the distribution the 
mode tending slightly greater than the mean. However, will seen 
that the normal curve which has been fitted close approximation the 
observed distribution. Referring Table summary K-values can 
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Table and Fig. show remarkable result—that the values are 
nearly the same regardless the class phenomenon from which they are 
determined. necessary examine them little further see the differ- 
ences are significant. will noticed that the standard deviation from 
all the values little larger than 
that for most the values for the 
groups. When the standard devia- 


tion the means groups cal- 
times the standard deviation the 
values This suggests that 
does vary slightly with different 
phenomena. mean value 


0.81 from tree rings, and 
this derived from four sets trees, 
each set being from one four dif- 
are 0.84, 0.80, 0.81, and 0.81, 

that would appear that the de- 
parture the trees from the mean 

all the phenomena significant. Similarly, the difference between 
the varves Lake Saki (the Crimea) and those Tamiskaming (Ontario, 
Canada) and Moen (Norway) greater than would expected chance 
distribution. Tamiskaming and Moen are significantly different. The rainfall 


TABLE K-VALUES 


Phenomena No. Mean 

Varves (Tamiskaming, Ont., Canada, and Moen, Norway)......... 0.77 0.094 
Sunspot numbers and wheat prices (combined miscellaneous 


observations have been examined detail, with the result that the means for 
fifteen stations with some values each station not differ greatly 
from random distribution. 

The general conclusion that, long series several hundred years 
available for phenomenon, probably best take the mean value 
from these indication for the future. However, only short period 
from 100 years 200 years less available, the best plan take the mean 
value for from all available material. The mean adopted for this purpose 
0.72 instead 0.73, order give little more weight the more precise 
measurements relating rivers, rainfall, and temperature. The mean value 
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that the mean value quantity that has ranged from 0.46 0.96, and, 
therefore, some margin safety must allowed. 

The probable error single value R/o can derived differentiating 
with respect and using the data Table which yield probable 
error 24%. Confirmation this value obtained taking the values 
for rainfall, river discharge, temperature, and pressure for values 
between years and 120 years, and correcting them period 100 years. 
The probable error calculated from these corrected values Fig. 10, 
introduced subsequently, shows how can reduced 23%, reducing 
the 0.04 

Fig. the mean representing all the lines represented 


0.72 


which means that R/o increases more rapidly with the case natural 
phenomena than the case chance events for which the power 
The value for 100 from Eq. 16.8. the computations previ- 
ously cited? the equation 1.65 was used, which gives R/o 16.5 
100. 

phenomenon exactly periodic, with constant amplitude, will have 
definite maximum which will reached the end each complete period. 
this case will not vary the manner already found for natural phenomena. 
The phenomena considered this paper, however, are not this type. 
the case sunspot numbers (which are the nearest approach) the periodic part 
does not have constant amplitude, and has irregular variations superposed. 
These irregularities cause increase increases and the mean value 
0.75, with indication that the phenomenon different from the many 
others considered. 

The very striking result from the investigation the small variation the 
index over such wide range phenomena, and hence the possibility its 
utility many fields not mentioned this paper. 


VARIATION THE MEAN AND THE STANDARD DEVIATIONS 
DERIVED FROM SHORT-TERM RECORDS 


When records natural phenomenon extend over long periods there are 
considerable variations both means and standard deviations from one period 
another. The primary practical interest this paper water supply 
and, hence, rainfall. Unfortunately 200 years about the longest series 
observations obtainable for rainfall, and few stations reach 150 years. Some 
these records are given Table They are divided into periods 
years years illustrate the recorded variations. This division chosen 
because quite frequently water supply hydroelectric projects have 
based such short records since nothing more available. The reservoirs 
the Great African Lakes are case point. 

The average range variation the mean values the rainfall Table 
about 14%, whereas the standard deviations vary the average 23%. 


st 
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This more than would expected from the ordinary theory for random 
distributions. The same true the case temperature. 

The uniformity the values shown the preceding section encourages 
the examination some the long-period records, summary which 


TABLE MEANS AND STANDARD DEVIATIONS 
THE CASE RAINFALL 


STANDARD DEVIATION 
(years) (years) Maximum Minimum Maximum Minimum 
Philadelphia, Pa........ 126 42.7 41.4 7.0 5.6 
Boston, Mass........... 128 38.2 5.0 
Madras, India.......... 133 48.4 15.6 14.3 
151 886 777 17.1 15.7 
Stockholm, Sweden..... 161 57.8 38.3 13.5 7.8 
Milan, Italy............ 171 104 21.2 15.1 
Zwanenbourg, the 
Netherlands.......... 211 13.1 10.8 


compiled Table 10. Table 10(a) the variability means for years 
and 100 years compared with what would expected the frequency were 
random, and not affected any tendency for high low values grouped 
together. 

The computed value the standard deviation mean for years 
(Col. obtained taking the standard deviation obtained from all the 


TABLE MEANS AND STANDARD DEVIATIONS 


MEANS DEVIATIONS 
Phenomenon years, MEANS PERIODS 
Varves: 
Tamiskaming, Ont., 
Ring Thickness: 
Roda gage (free 


individual observations for the period given Col. and dividing 50. 
Cols. Table 10, are the ratios Cols. 3/Cols. appears from Table 10(a) 
that the means are much more variable than they would the case 
random distribution, the standard deviation 50-year mean being 2.5 times 
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great would for random distribution. Similarly, the standard 
deviation 100-year mean 3.2 times great random distribution. 
This implies that, over long period, the maximum value likely higher 
and the minimum lower than would predicted from the application the 
ordinary theory probability the records from short period. 

Table shows the variability the standard deviation for periods 
years and 100 years those cases where the records cover long time. 
distribution normal and randomly distributed chance events, the 
standard deviation the standard deviation This the value 
called (Cols. actual value (Cols. found the usual 
way from the standard deviations for all the 50-year periods, and Cols. 
contain the ratios, Cols. 3/Cols. Again appears that the standard 
tion for these natural phenomena more variable than would the case 
random distribution. 

many natural phenomena have nearly normal fre- 
quency distribution, this only the case when their order occurrence not 
considered. Consequently, their representation random distribution 
the usual way only first approximation. The tendency occur groups 
makes both the mean and the standard deviation computed from short periods 
years more variable than the case random distributions. Thus the odds 
against the occurrence some particular extreme value (say, high flood) are 
likely less than those calculated from short period records normal 
probability theory. 

was thought that might possible derive expression for the 
effect grouping, involving the correlation coefficient between successive 
values the variable. Experiments with the natural number ---, and 
—1, —N, however, showed that this was not the case. The greatest 
value obtained placing the positive numbers one group and the 
negative another and independent the the groups. The correla- 
tion between successive numbers depends the order within the groups. 
the foregoing example, the coefficient correlation may vary from while 
remains constant N(N The numbers can also arranged give 
small values with correlation coefficients varying from —1. 


10. GUARANTEE Less THAN 
THE MEAN DISCHARGE 


the discharge from reservoir over period kept steady the mean, 
the range the curve progressive sums departures from the mean taken 
the storage required keep this constant discharge. The final sum the 
departures zero. departures are taken from base less than the mean, this 
corresponds the case which the reservoir discharge not allowed fall 
below this base value. 

The storage required this case the greatest accumulated deficit. This 
illustrated Fig. Lake Albert discharges, which the summation 
curve, curve accumulated departures from base 24, plotted, axis 
Fig. and the summation curve for departures from the mean referred 
axis OX. 
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departures are computed from new base below the mean, 
explained previously, the summation curve must referred the inclined 
that, the inclination the axis changes, different points the summa- 
tion curve may become maxima, minima, and there simple connection 
between the storages required guarantee the mean and guarantee some 
lesser discharge. This clearly shown Fig. which the result tossing 


Number Throws 


coins 100 times shown the form curve the continued sums the 
number heads minus the number tails. 

the end the 100 tosses, heads exceeded tails twelve, that the 
departures heads minus tails are reckoned from base 0.12 below the mean. 
This corrected Fig. drawing the axis which the base for the 
summation curve departures heads minus tails from their mean. Simi- 
larly, the axes OX, and can drawn for the summation curves de- 
partures from the bases 0.1 and below the mean. The two points 
are the maximum and minimum the summation curve heads minus tails 
and also the summation curve departures heads minus tails from its 
mean. The range the first curve greater than the range the second. 


Mean 
thmetic 
from 
Mea, 
-40 
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The maximum accumulated deficit with the mean base given the points 
ab, but when base 0.1 selected below the mean the maximum accumulated 
deficit given for base 0.2 below the mean, The storage 


Media 
required maintain draft below the mean can computed arithmetically; 
but extremely laborious and far the quicker and more reliable method 
the graphic one. 
VALUES 
Magnitude 
Rivers, RAINFALL, AND TEMPERATURE 
0.48 0.35 0.27 
Log —0.32 —0.46 —0.57 


0-0 
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11. BETWEEN AND STORAGE 
FoR CHANCE EVENTS 


Let the storage required maintain draft less than the mean (M) 
Values (the maximum deficit) have been computed from the experiments 
with coins, probability cards, and drawing bonds, described. They have 
been calculated for every set 100 experiments and for values 
The ratio S/R has been tabulated and mean values are given 
Table 11. 

The results are also given Fig. 7(a) which log S/R plotted against 
the relation being 


little better fit obtained from the equation— 


which shown Fig. Neither formula should used outside the range 
values from which was computed. 

Considerable research was devoted finding theoretical relation between 
and but without success. extension the method Section will 
enable the mean range the summation curve departures from base below 
the mean found, but does not result This mean range 
not what wanted, however; the maximum accumulated deficit the 
entire period, any part that wanted. The latter complicated 
function gains minus losses, for which algebraic expression has been found. 

interesting point the case coins, cards, etc., that, the mean range 
called when departures are measured from base less than the mean, 
then the ratio approximately equal S/R. 


Maximum Dericit (MEAN VALUES) 


Magnitude 

0.5 0.52 0.53 0.55 0.6 0.63 0.66 0.68 0.7 0.88 

0.34 0.32 0.24 0.24 0.23 0.20 0.20 0.09 S/R 

0.21 0.15 0.12 S/R 
—0.68 —0.82 —0.92 log (S/R) 
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NaTURAL PHENOMENA 
The results number computations the storage (S) required 
guarantee minimum draft (B) less than the mean have been compiled. 
Table gives the mean values these detailed computations. Curves 


River 
Temperature 


Rainfall 
Temperature 


River RAINFALL, AND TEMPERATURE DETERMINE 


the exponential and square root forms have been fitted these mean values 
and the results are shown Figs. and far closeness fit con- 
cerned, over the range observations, there significant difference be- 
tween one type curve and the other. some future time, may perhaps 


0-0 


S/p 


possible decide that one the other has some theoretical justification, 
but this has not far been possible. The results the curve fitting may 


studied Table 13. The curves that have been fitted the means all the 
observations are shown Fig. 10. 


N O 
(a) 
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The equations show very little variation among the three groups phe- 
nomena. For example, 0.25, S/R from the exponential curve 
for each group 0.47, 0.47, and 0.42 and from the square root curve for each 
group 0.46, 0.49, and 0.42. These results, again, illustrate striking 
uniformity. 


TABLE 13.—Comparison SHOWING RELATION BETWEEN 


Fig. No. Media Log (S/R) S/R 


Considering the curves based the means all the observations, which 
are given Fig. 10, will seen that spite the difference form, the 
foregoing expressions give practically identical curves over nearly all the range. 

The important result shown 
maximum value (the mean) makes 
much greater proportional reduc- 
tion the storage required main- 
tain it. For example, reduction 
0.1 diminishes the storage 
about 35%. This fact provides 
very practical method applying 
factor safety compensate 
the natural variations 


FROM DIFFERENT SOURCES 
ALL OBSERVATIONS ANALYZED 


with the storage one variable 
supply, for example, the water flowing out Lake Albert. was shown 
(see Eq. 13) that, the average, the capacity required produce steady 
discharge equal the mean over period years, stream whose 


0.72 
two streams are involved, with variations whose discharges have normal 
distribution, then and are their standard deviations, and the 


standard deviation their combined flow, can shown that 


which the coefficient correlation between the two discharges. 
follows from Eq. that, for given value 


= 


RESERVOIR CAPACITY 


Except this case the necessary capacity equalize the discharge the main 
stream less when all one reservoir the main stream than when 
reservoirs each the tributaries. possible under some circum- 
stances for the discharge the main stream equalized reservoir 
tributary. 

14. APPLICATIONS 

mentioned Section the investigation arose out the need control 
the flow the White Nile reservoirs Lake Albert and Lake Victoria. 
The scheme favored Egypt was construct single large reservoir Lake 
Albert such capacity equalize the outflow over century. Such 
scheme requires the simplest type regulation, but was not acceptable 
Uganda because the large area land which would flooded for long 
periods. The Uganda Government proposed the combination large storage 
reservoir Lake Victoria with small regulating reservoir Lake Albert. 
Later regulator the exit Lake Kioga Uganda was seen necessary 
part this combination eliminate the very variable effect Lake Kioga 
the Victoria Nile. The regulation the combination more complicated and 
has been the subject much study. 

the more simple plan for Lake Albert the capacity the reservoir would 
large compared with the natural fluctuations the lake. The method 
working have been build the reservoir years years before 
was required complete action, and, during this time, fill partly full, 
leaving sufficient reserve capacity unusual floods that might occur 
after the project was working full capacity and which would then damage 
the Sudan. The simplest working the reservoir would send down 
constant annual discharge slightly less than the mean, distributing this dis- 
charge required through the year. This operation would lead occasional 
years which discharge more than the average would necessary 
avoid overfilling the reservoir. 

This procedure was adopted the report? for making estimates the value 
the projects and the size the works required. these estimates the 
guaranteed discharge was less than the mean discharge twice the probable 
error the mean, plus further small reduction. The actual procedure would 
depend circumstances, beginning, doubt, with what was thought 
safe draft. the reservoirs showed signs reaching high levels, the draft 
would increased and, the contents diminished considerably, the draft 
would reduced. During the early years regulation, before the full draft 
was required, experience would gained which would enable sound scheme 
regulation drawn based more information than now available. 

modification this simple regulation would include what the writer has 
called which excess the Ethiopian tributaries can 
virtually stored reducing the outflow from Lake Albert. 

Lake Victoria, the main storage reservoir, has the advantage that its area 
12.5 times that Lake Albert. Therefore, much lower dam will give the 
required storage; but the regulation complicated the fact that the natural 
range rise and fall the lake greater than the extra range required for 
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regulation. This would not matter there were plenty room for regulation, 
but the interests the inhabitants around the lake place definite limit 
the height which the lake can raised, and the natural range the lake 
plus the extra effect regulation must fitted within the allowable limits. 
This requirement fixes the capacity needed Lake Albert, where extra storage 
over that required for the annual regulation necessary provide against 
unusual floods. The regulation the two lakes combination, therefore, 
departs from the simple one steady annual discharge out Lake Victoria 
and steady annual discharge out Lake Albert, which would almost auto- 
matic. The more complicated regulation will depend for success upon close 
study rainfall, runoff, discharges tributaries, and trends lake levels. 
However, the time that will elapse before the project working will provide 
additional data, and experience regulation. 

fortunate that the Nile supply comes from two sources—Ethiopia and 
the Great African Lakes—and that, far available information goes, the 
correlation between their contributions small. unlikely, therefore, that 
both sources would suffer extended periods drought together, and this can 
provide factor safety against unusual sequence low years the Lake 
Plateau which includes Uganda and parts Kenya, Tanganyika, and the 
Belgian Congo. take advantage this the scheme projects reported? 
must considered whole and worked combination. Should 
unusual succession low years the Lake Plateau cause serious decline 
the stored water, would necessary make full use virtual storage. 
Then, beginning store little earlier the main Nile reservoirs, the quota 
from the lakes could reduced and the possibility dangerous depletion 
the storage would removed. 


15. 


This paper has been mainly confined theoretical matters, and the prac- 
tical outcome therefrom has been reported 
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APPENDIX. NOTATION 


The following letter symbols have been adopted for use this paper and 
the discussion: 


annual draft from storage stream, less than the mean annual 
flow the stream; 
number combinations different things taken time; 
departure: 
average difference between the range the continued sums 
continued sums, cumulative totals; 
base natural logarithms; 
number gains; 
total number instances which the losses exceed the gains; 
log N/2’ 
number losses; 
mean annual draft water from storage; 
number coins any toss (see C); 
number annual discharges recorded; period observation; also, 
number tosses, set coins; number departures from 
mean; 
Nm; 


index 


discharge, annual; 


range, from maximum minimum, the curve obtained plotting 
the cumulative totals the departures from the mean 
the annual discharges stream, taken order; 

number heads tosses coin; also; number heads; 
(see C); coefficient correlation between two variables; 

storage required maintain draft less than the mean maxi- 
mum storage deficit; 

abscissas; 

ordinates; 

“sum and 

standard deviation: 
the number heads minus the number tails—that 

is, 2m. 


y= 


q 

j 

] 
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